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ABSTRACT 

We develop a density tensor hierarchy for open system dynamics, that 
recovers information about fluctuations (or "noise") lost in passing to the re- 
duced density matrix. For the case of fluctuations arising from a classical 
probability distribution, the hierarchy is formed from expectations of prod- 
ucts of pure state density matrix elements, and can be compactly summarized 
by a simple generating function. For the case of quantum fluctuations arising 
when a quantum system interacts with a quantum environment in an overall 
pure state, the corresponding hierarchy is defined as the environmental trace 
of products of system matrix elements of the full density matrix. Whereas 
all members of the classical noise hierarchy are system observables, only the 
lowest member of the quantum noise hierarchy is directly experimentally 
measurable. The unit trace and idempotence properties of the pure state 
density matrix imply descent relations for the tensor hierarchies, that relate 
the order n tensor, under contraction of appropriate pairs of tensor indices, 
to the order n — 1 tensor. As examples to illustrate the classical probability 
distribution formalism, we consider a spatially isotropic ensemble of spin- 
1/2 pure states, a quantum system evolving by an Ito stochastic Schrodinger 
equation, and a quantum system evolving by a jump process Schrodinger 
equation. As examples to illustrate the corresponding trace formalism in the 
quantum fiuctuation case, we consider the tensor hierarchies for coUisional 
Brownian motion of an infinite mass Brownian particle, and for the the weak 
coupling Born-Markov master equation. In different specializations, the lat- 
ter gives the hierarchies generalizing the quantum optical master equation 
and the Caldeira-Leggett master equation. As a further application of the 
density tensor, we contrast stochastic Schrodinger equations that reduce and 
that do not reduce the state vector, and discuss why a quantum system cou- 
pled to a quantum environment behaves like the latter. The descent relations 
for our various examples are checked in a series of Appendices. 
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1. Introduction 



Increasing attention is being paid to the dynamics of open quantum systems, that is, 
to quantum systems acted on by an environment. Such systems are of interest for studies of 
dissipative phenomena, decoherence, backgrounds to quantum computers and to precision 
measurements, and theories of quantum measurement. A principal tool in studying open 
quantum systems is the reduced density matrix, obtained from the pure state density matrix 
by tracing over environment degrees of freedom, or in stochastic models where the environ- 
ment is represented by a noise term in the Schrodinger equation, by averaging over the noise. 
As is well-known, this transition from the pure state density matrix to the reduced density 
matrix is not one-to-one, since information about the total system is lost. For example, in 
stochastic models, there is known to be a continuum of different unravelings, or pure state 
density matrix stochastic evolutions, that yield the same master equation for the reduced 
density matrix. The question that we investigate here is the extent to which one can form 
objects that refer only to the basis vectors of the system Hilbert space, but that nonetheless 
recapture information that is lost in passing to the reduced density matrix. In the first 
part of this paper (Sections 2 through 5), we discuss classical noise arising from fluctuations 
defined by classical probability distributions. In the second part (Sections 6 through 9), we 
give an analogous discussion of quantum noise, which appears in the physically important 
case of a quantum system coupled to a quantum environment in an overall pure state. We 
also give an extension, making contact with the discussion of the first part, to the case in 
which the overall system is in a mixed state superposition of pure states. The final section 
contains a discussion of quantum measurements that relates the material in the first and 
second parts. 
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For the case of classical probability distributions, a relevant discussion appears in 
Chapter 5 of the book The Theory of Open Quantum Systems by Breuer and Petruccione 
[1], following up on earher papers by those authors [2], by Wiseman [3] and by M0lmer, 
Castin, and Dalibard [4]. In simplified form, Breuer and Petruccione introduce an ensemble 
of pure state vectors \ipa), each drawn from the same system Hilbert space Hs, with each 
vector assumed to occur in the ensemble with probability Wa, S^Wa — 1. Measurement 
of a general self-adjoint operator R for a system prepared in \ipa) typically gives a range of 
values, the mean of which given by {ipalRli^a) ■ The mean or expectation over the ensemble 
of pure state vectors is then given by 

a 

with p the mixed state or reduced density matrix defined by 

a 

Breuer and Petruccione point out that there are three variances that are relevant. The 

variance of measurements of R over all pure states in the ensemble is given by 

Var(i?) = Trp{R - TrpRf = Trpi?^ _ (Trpi?)^ . (2a) 

This can be written as the sum of two non-negative terms, 

Var(i?) = Vari(i?) + Var2(i?) , (26) 

with Vari(it!) the ensemble average of the variances of R within each pure state of the 
ensemble, 

ys.Tl{R)=J2M{^a\R'm-{i^a\R\i^af] , (2c) 
a 
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and with Var2(it!) the variance of the pure state means of it! over the ensemble, 

Var2(i?) = J]^/;a(^a|i?|V'a)'-E^a(V'a|i?|V'a)]' ■ (2rf) 

a a 

Thus, Vari(i?) is an ensemble average of the quantum variances of while Var2(i?) is a 
measure of the spread of the average values of it! resulting from the statistical properties of 
the ensemble. As Breuer and Petruccione note, neither of the subsidiary variances Vari 2 
can be expressed as the density matrix expectation of some self-adjoint operator. 

Our aim in the first part of this paper is to extend the formalism of ref [1] by utilizing 
a density tensor hierarchy, which captures the statistical information that is lost in forming 
the reduced density matrix of Eq. (lb). A density tensor, defined as an ensemble average 
of density matrices, was first introduced by Mielnik [5], and was applied to discussions of 
density functions on the space of quantum states and their application to thermalization of 
quantum systems by Brody and Hughston [6]. These papers, in addition to introducing the 
concept of a density tensor which is developed further here, also contain the important result 
that in the case of a continuum probability distribution, the density tensor hierarchy gives 
all of the information needed to reconstruct the probability function Wa- In particular, the 
variances Vari 2 for any observable, and more general statistical properties of the ensemble 
as well, can be expressed as contractions of density tensor matrix elements with appropriate 
matrix elements of the observable (s) of interest. 

The basic construction of the density tensor hierarchy corresponding to a classical 
probability distribution {wa} is given in Sec. 2. Here we generalize the reduced density 
matrix of Eq. (lb) to a density tensor, formed by taking a product of pure state density matrix 
elements, and averaging over the ensemble of pure states. When the ipa are independent of 
a, this tensor reduces to an n-fold product of reduced density matrices, and so the difference 
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between the density tensor and this product is a measure of the statistical fluctuations in the 
ensemble. In the generic case of non-trivial dependence of ipa on a, there are some general 
statements that can be made. First of all, the order n density tensor is a symmetric tensor 
in its pair indices, and it can be considered as a matrix operator acting on the n-fold tensor 
product of the system Hilbert space Hs with itself. The symmetry of the density tensor 
allows construction of a generating function that on expansion gives the density tensors of 
all orders. Additionally, as a consequence of the unit trace and idempotence conditions 
obeyed by the pure state density matrix, the density tensor hierarchy satisfles a system of 
descent equations, relating the order n tensor to the order n — 1 tensor when any row index is 
contracted with any column index. We show that the variances Vari 2 deflned by Breuer and 
Petruccione can be expressed in terms of appropriate contractions of density tensor elements 
with operator matrix elements. 

In subsequent sections we develop some concrete applications of the general formal- 
ism for classical probability distributions. In Sec. 3, we consider an isotropic ensemble of 
spin-1/2 pure state density matrices, construct the density tensors through order 3, ver- 
ify the descent equations, and calculate the generating function. In Sec. 4 we apply the 
formalism to a quantum system evolving under the influence of noise as described by a 
stochastic Schrodinger equation, with the ensemble deflned as the set of all histories of an 
initial quantum state under the influence of the noise. Assuming white noise described by 
the Ito calculus, we give the dynamics of the general density tensor in terms of the general 
unraveling of the Lindblad equation constructed by Wiseman and Diosi [7], and show that 
the order two and higher density tensors distinguish between inequivalent unravelings that 
give the same reduced density matrix (i.e., the same order one density tensor). In Sec. 5 
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we develop an analogous formalism for the case of jump (piecewise deterministic process) 
unravelings of the Lindblad equation. 

We turn next to an analysis of a quantum system coupled to a quantum environment, 
rather than to an external classical noise source. Here, one is confronted with the problem of 
discussing the system dissipation associated with the system-environment interaction within 
a single overall pure state of system plus environment (or in a thermal state that is a 
weighted average of such pure states). Typically, in master equation derivations, the system- 
environment interaction^ H has vanishing expectation in the environment, but its square 
does not have a vanishing expectation, because the environment is not in an eigenstate of 
H. The associated variance is then a measure of quantum fluctuations associated with the 
environment state, and is the source of quantum "noise" driving the system dissipation. 
Our aim in the second part of this paper is to generalize the formalism of the first part 
to recapture information about this noise that is lost in the passage to the system reduced 
density matrix. We do this in Sec. 6 by defining a density tensor hierarchy as the trace over 
the environment of a product of environment operators constructed as the system matrix 
elements of the total density matrix. Unlike the classical noise construction, which uses only 
the system density matrix, the construction in the quantum noise case requires knowledge of 
the full system plus environment density matrix, and so (except for the order one case) does 
not give a system observable. It is nonetheless computable in any theory of the system plus 
environment, and is of theoretical, rather than empirical, interest. Because the environment 

^ What we call H is usually denoted by Hj in the open systems hterature. To avoid 
confusion, all other Hamiltonians will carry subscripts, e.g., Hs and Hs for the system and 
environment Hamiltonians, i^TOT for the total Hamiltonian, etc. 
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operators entering the construction are non-commutative, this hierarchy is no longer totally 
symmetric in its system index pairs, but by the cyclic permutation property of the trace, it 
is symmetric under cyclic permutation of the system index pairs. Also, because the system 
trace of these environment operators gives only the reduced environment density matrix, 
rather than unity, there is in general no descent equation associated with taking this trace. 
However, when indices of adjacent system operators are contracted, one gets the square of 
the overall density matrix, and so there remains a set of descent relations connecting the 
order (n) tensor to the order {n — 1) tensor. Finally, in the case of thermal (or other mixed) 
overall states, we define the appropriate tensor as a weighted sum of pure state tensors, in 
analogy with the definition of Sec. 2. 

In subsequent sections, we give applications of the trace hierarchy formalism to 
several classic problems discussed in the theory of quantum master equations. In Sec. 7 we 
consider the quantum Brownian motion (and resulting decoherence) of a massive Brownian 
particle in interaction with an independent particle bath of scatterers. In Sec. 8 we discuss 
the tensor hierarchy corresponding to the weak coupling Born-Markov master equation, 
and it specialization to the quantum optical master equation. Finally in Sec. 9, we give 
an analogous discussion for the Caldeira-Leggett model of a particle in interaction with a 
system of environmental oscillators. 

We conclude with a discussion that bridges the considerations of the classical noise 
and the quantum noise cases. In Sec. 10, we contrast two different Ito stochastic Schrodinger 
equations, both of which have the same Lindblad, but only one of which leads to state vector 
reduction. We relate this to the fact that the equation giving the time derivative of the 
stochastic expectation of operator variances involves the order two density tensor, which 
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differs for the two cases. We discuss the analogous equation for the time dependence of the 
variance of a "pointer operator" in the case of a quantum system coupled to a quantum 
environment, and show why this does not lead to state vector reduction. Thus we see no 
mechanism for quantum "noise" in a closed quantum system plus environment to provide a 
resolution of the quantum measurement problem. 

2. The density tensor for classical noise and its kinematical properties 

We proceed to establish our notation and to define the density tensor hierarchy in 
the classical noise case. We denote the pure state density matrix formed from the unit 
normalized state |'0ck) by Pq,, with 

Pa = |^a)(^a| , (3a) 

and its general matrix element between states \i) and \j) of Tis by 

Pa;ij = (ilpalj) ■ (36) 

The unit trace condition on states that 

Tipa = (^/'al^a) = 1 , (3c) 

and the idempotence condition on states that 

pI = \'ipa){lpa\\'^a){lpa\ = \lpa){lpa\ = Pa ■ (3rf) 

We now define the order n density tensor by 

Pii]i,i2i2 inin ~ ^ '^aPa;hji Pa;hj2 ■ ■ ■ Pa;injn — -^[Pa;iijiPa;i2j2---Pa;injn] ) (^o) 



with E[Fa\ a shorthand for 

a 

Since 

a a 

we see that this is just the \i) to \ j) matrix element of the reduced density matrix p defined 
in Eq. (lb), 

and so the density tensor of Eq. (4a) is a natural generalization of the usual reduced density 
matrix. When the states \ipa) are independent of the label a, the definition of Eq. (4a) 
simplifies to 

(n) (K \ 

and so the difference between Eq. (4a) and a product of reduced density matrix elements 
is a reflection of the statistical structure of the ensemble. Since the factors within the 
expectation E[..] on the right of Eq. (4a) are just ordinary complex numbers, the density 
tensor is symmetric under interchange of any index pair iiji with any other index pair imjm- 
Consequently, we can define a generating function for the density tensor by 

oo 
n=0 

where repeated indices i,j are summed. It will often be convenient to abbreviate Pa-ijdij by 
Pa • a, so that the generating function becomes in this notation G[a\ — £J[e^"'"]. 

Although the density tensor for n > 1 is not an operator on Tis, it clearly has the 
structure of an operator on the n-fold tensor product Hs <E) Hs ■■■ ® Hs- Motivated by 
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this, we will often find it convenient to write the definition of Eq. (4a) as 



= E 



n 

.C=l 



Pa-/ 



(5e) 



with each factor pa-/ acting on a distinct factor Hilbert space Hs-e in the tensor product 
nr=i '^S;i- One can pass easily back and forth from this notation to one in which the system 
matrix indices are displayed explicitly. 

Let us consider next the result of contracting any row index ii with any column index 
jk- There are two basic cases: (i) one can contract a row index ii with its corresponding 
column index ji, and (ii) one can contract a row index ii with a column index jk with k ^ I. 
Since the density tensor is symmetric in its index pairs, it suffices to consider only one 
example of each case, since all others can be obtained by permutation. For the contraction 
of ii with ji we find 

^iijiPiiji,i2h,-,injn ~ ^[i'^P)Pa;i2j2---Pa;injn\ — ^[Pa;i2j2---Pa;injn] — Pi2j2j..,injn ' (^^) 

where we have used the unit trace condition of Eq. (3c). For the contraction of ji with ^2, 
we find 



^jii2Phji,i2j2,--;injn "^[(P ) a;iij2 ■ ■ ■ Pa;i„j„] E[Pa;hj2Pa;i3j3---Pa;injn] Piij2,i3j3,-,injn ' (^^) 

where now wc have used the idempotence condition of Eq. (3d). As an illustration of how 
this works when all possible index pair contractions are considered, we give the complete set 
of contractions reducing the second order density tensor to a first order density tensor. 



.(2) 



(7a) 

"jl«2FllJl,J2J2 f^llJ2 ' 

5- ■ p^^- ■ ■ ^p^^- 

"J2J1PJ1J1,12J2 ^«2J1 
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Referring to the generating function of Eq. (5d), the general descent equations can be sum- 
marized compactly by the two identities, 

Smr^§^ =^[(Trp„)e^--"-] = G[a,,] , 

OOj-priT (7&) 

When the density matrix p used to define the density tensor is a mixed state density matrix, 
the trace descent relation of Eq. (6a) is unchanged, while the indempotency relation of 
Eq. (6b) relates the contraction an order (n) tensor to an order {n — 1) tensor in which one 
factor p is replaced by p^; this is not a member of the original hierarchy, but still gives a 
useful relation for checking calculations. 

To conclude this section, let us return to the variances introduced by Breuer and 
Petruccione. In terms of the order one and order two density tensors, we evidently have 

Vari(i?) ^p\l]^{R'^)j^i^ - pf^j^,i^j^RhnRj2i2 , 

Var2(i?) =f^,,,,,Rj,nRj,i, - {pfXRn^.f , 

Var(i?)=plJ]^(i?^),,,-(p«%,)^ , 
with Rji — {j\R\i). Clearly, other statistical properties of the ensemble are readily expressed 
in terms of the density tensor hierarchy. For example, the ensemble average of the product 
of the expectations of two different operators R and 5" is given by 

^W„(V'a|^|V'a)(V'a|'5|V'a) , (86) 

a 

which can be used, together with information obtained from p^^^ , to calculate the covariance 
and correlation of R and 5*. 
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3. Isotropic spin-1/2 ensemble 



As a simple example of the density tensor formalism, let us follow Breuer and Petruc- 
cione [1] and consider the case of an isotropic spin-1/2 ensemble. Let v be a vector in three 
dimensions, and consider the ensemble of spin-1/2 pure state density matrices 



Pi.^-\i^ + v-a) , (9a) 



with a — ((7^, (7^, (7^) the standard Pauli matrices, and with a uniform probability distribution 
of V over the unit sphere 1^ | = 1 specified by 

w{v)^^5{\v\-l) . (96) 

(Clearly, v has the same significance as the label a used in the preceding section.) Defining 

E[P{v)] = J Svw{v)P{v) , (10a) 

a standard calculation gives 



E[l\ = 1 , E[vsVt\ = ^5st , ... , (106) 



with all averages of odd powers of v vanishing. From Eq. (9a), we have 



and the general density tensor over this ensemble is defined by 
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Prom Eq. (10b), the first three tensors in this hierarchy are now easily found to be 



(2) U . 



(3) 1 
■ . g 



^iiji^i2j2^i3j3 ~^ 2 {^hji'^i2j2 ' ^hh ~^ ^i2j2^hji ' ~^ ^hh^hji ' ^^232) 



(12) 

Using the relations Tra — and (?^)y = 35^, it is now easy to verify that the descent 
relations of Eqs. (6a) and (6b) are satisfied by Eq. (12). 

For the isotropic spin-1/2 ensemble, the generating function of Eq. (5d) becomes 

G[aij] = Ele"^^'^'"'^] , (13) 

with p{v )ij given by Eq. (11a). Defining the vector A by 

A= ^aijUij , (14a) 

a simple calculation gives 

1 Ar^ ,sinhU| ,1^ (A'^f , , 

GM=exp(-IVa)^^ = exp(-IVa)[l + — + ^ + ...] , (146) 

from which one can read off the values of the low order density tensors given in Eq. (12). 
The verification of the descent relations of Eq. (7b) for the generating function of Eq. (14b) 
is given in Appendix A. 

4. Ito stochastic Schrodinger equation 

We consider next a state vector 1-0) with a time evolution described by a stochas- 
tic Schrodinger equation, which is a frequently used model approximation to open system 
dynamics. In this case the state vector and the corresponding pure state density matrix 
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p — ^re implicit functions of the noise, which takes a different sequence of values for 

each history of the system. In the notation of Sec. 2, the different histories are labeled by 
the subscript a, and the expectation of Eq. (4b) is an average over all possible histories. It 
is customary, however, in discussing stochastic Schrodinger equations to omit the subscript 
a, treating the history dependence of p as understood. So in this context, the definition of 
Eq. (4a) becomes 

with E[..] the stochastic expectation, and the generating function G[aij\ takes the same form 
as given in Eq. (5d) but with the subscript a omitted. 

Our aim in this section is to derive an equation of motion for the generating function, 
which on expansion yields equations of motion for all density tensors taking as input 
the general pure state density matrix evolution constructed by Wiseman and Diosi [7] , that 
corresponds to a given Lindblad form [8,9] for the time evolution of the reduced density 
matrix p'^^^ — E[p\. We begin by recapitulating the results of ref [7]. The most general 
evolution of a density matrix p that preserves Trp = 1 and obeys the complete positivity 
condition is the Lindblad form 

dp — dtCp , (16a) 

with 

Cp = -i[HTOT, p] + Ckpcl - ^{4cfe, p} , (166) 

with {, } denoting the anticommutator, and with the repeated index k summed. The set of 
Lindblad operators Ck describes the effects on the system of the reservoir or environment that 
is modeled by an external classical noise. Wiseman and Diosi show that the most general 
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evolution of the pure state density matrix p for which E[dp\ reduces to Eqs. (16a) and (16b) 
takes the form 

dp = dtCp+\d(t))l(ip\ + \^){d(t)\ . (17a) 
Here \d(j)) is a state vector that is a pure noise term, so that 

E[\d(t>)]^Q , (176) 

that is orthogonal to |^), so that 

m<P)=0 , (17c) 

and that obeys 

\d(l)){d(l)\ = dtW . {17d) 

The operator W is the Diosi transition rate operator [5] given by 

W^Cp-{p,Cp} + pTr{pCp) 

(18) 

=(cfc - {ck))p{ck - {ck)y 

where (cfe) is a shorthand for the quantum state expectation {ip\ck\ip) = Trpc^. Although 
\d(l)){d(l)\ is completely fixed, Wiseman and Diosi show that \d(j))\d(/)) is free, with different 
choices for this and different phase choices for the corresponding to different pure state 
evolutions (or "unravelings" ) that yield the same evolution of Eqs. (16a) and (16b) for the 
reduced density matrix p. 

Wiseman and Diosi further show that \d(j)) can be parameterized by complex Wiener 
processes by writing 

\d<f>)^{ck-{ck)mdCk , (19a) 
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with 

E[d^,] = E[dC,] = (196) 

and with 

dCj{t)dCk{t) ^dtSjk 

(19c) 

d^j{t)d^k{t) =dtUjk , 

where Ukj — Ujk is a set of arbitrary complex numbers subject to the condition that the 
norm of the complex matrix u = [ujk] be less than or equal to 1. (See Eqs. (4.10) and (4.11) 
of ref. [7].) In terms of this parameterization of \d(f)), the pure state evolution of Eq. (17a) 
takes the form 

dp = dtCp + (cfc - {ck))pdCk + p(cfc - {ck)yd^k , (19rf) 

and the corresponding stochastic Schrodinger equation for the wave function is [7] 

d\il;) = - iH^dt\iP) + (cfe - {ckmM , 

1 , (19e) 
-iH^ = - iHTOT - ^{clck - 2{ck)*Ck + (cfc)*(cfc)) . 

We proceed now to use pure state evolution of Eq. (19d) to calculate the evolution 

equation for the generating function 

G[aij] = E[exp{pijaij)] . (20a) 

To calculate the differential of Eq. (20a), we use the Ito stochastic calculus rule for the 
differential of a function f{w) of a stochastic variable w, 

df{w) = dwf{w) + ]^{dwff"{w) . (206) 

Applying this to Eq. (20a), we get 

dG[aij\ = E[{dpmrarnr + ^dpmramrdppgapg) exp{pijaij)] . (20c) 
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Substituting Eq. (19d) for dp, and using Eqs. (19a-c), together with the Ito calculus rule 

E[dwf{w)] = 0, we get 

dG[aij] = dtE[(amr{jO-p)mr + ^amrapgCmr,pq) ex.p{pijaij)] , (21a) 
with the coefficient of the quadratic term in aij given by 

Cmr,pq — Cpq^jju- — dpmrdppq 

= {m\{ck - {ck))p\r){p\p{ck - {ck))^\q) 

+(m|p(cfe - (cfc))t|r)(p|(cfc - {ck))p\q) (216) 

+(m|(cfe - (cfe))p|r)(p|(Q - {ce))p\q)ulf 

+ {m\p{ck - {ck)y\r){p\p{ci - {ci)y\q)uke 
This expression can be rearranged by using the identity, valid for general operators A,B, 
general states |r), |m), and general pure state (idempotent) density matrix p, 

pA\r){m\Bp^ p{m\BpA\r) , (22a) 

giving an alternative result for Cmr,pq 

Cmr,pq — ^^mqPpr ~l~ ^^prpmq 

+ [(Cfc - {Ck))p]mqU*u[{ct - {ci))p]pr (226) 
+ [p(Cfc - {Ck))'^]prUu[p{ce - (Q))^]mg , 

where we have used Eq. (18) defining the operator W , and where we use the subscript 
notation of Eq. (3b) for matrix elements, so that in general A^r = {m\A\r). 

Prom the evolution equation of Eqs. (21a,b) and (22b) for the generating function, 
by expansion in powers of a we can read off the evolution equation for the general density 
tensor of order n. Employing now the condensed notation of Eq. (5e), in which matrix 
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indices are not indicated explicitly, we have 

n 

dp^'^) =dtE[Y^{p,...p^),{Lp), 

(23a) 

+ ^ {Pl--Pn)tmCtm] ■ 

e<m=l 

Here {pi---Pn)e denotes the product YYj=i Pj with the factor p^ omitted, and similarly, {pi...pn)e.m 
denotes the product 11^=1 Pj "^i^^ factors p(, and pm omitted.^ The coefficient Ctm is 
given by 

Cim = Cmi =[(Cfe - {Ck))p]l[p{Ck - {Ck))'^]m 
+ [p{Ck - {Ck))%[{.Ck - {Ck))p]m 

(236) 

+ [(Cfc - (Cfc))p]4(Cfc - (Cfc))p]^M*^ 

+ [piCk - (Cfe))^]4p(Cfc - {c-k)y]rnUk-k , 

which corresponds in an obvious way to Eq. (21b) when matrix elements are written explicitly 
between states (m| and |r) in the Hilbert space labeled by i, and between states (p| and j^) 
in the Hilbert space labeled by m. (No relation is implied between the m used as a state 
label, and the m used as a Hilbert space label.) Since in Eq. (23a), which depends 
through the terms involving u^k on the choice of unraveling, is multiplied by two powers of 
a, it does not contribute to the evolution equation for the reduced density matrix p^^\ So 
as expected, the reduced density matrix evolution is given solely by the Lindblad term and 
is independent of the choice of unraveling. Higher density tensors p*^"), with n > 2, have 
evolution equations that receive contributions from Cem, and so contain information that 
distinguishes between different unravelings of the Lindblad evolution. 

As a simple illustration of how the tensors p^"^ for n > 2 distinguish between different 
unravelings, let us consider the case of real noise, d^k = d^k^ for which ujk = Sjk, and with 
^ For n = 1, (pi)i = 1 and {pi)em = 0, while for n = 2, (pip2)i2 = 1- 
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a single Lindblad ci, which we choose as either ci — A or ci — iA, with A a self-adjoint 
operator. Both choices of ci lead to the same Lindblad, since C is invariant under rephasing 
of Cfc, but through the Uj^^. terms they lead to different expressions for Cmr,pq- When = A, 
we find from Eq. (21b) 

C^n^.,p.=(m|{A- (A),p}|r)(p|{A- (A),p}|g) 

(24a) 

Hm\[pAP.M\\r){p\[pAP.MM , 
while when = lA, we have instead 

Cmr,p,^-{m\[A,p]\T){p\[A,p]\q) . (246) 

We will return to this example in Sec. 10. 

Using the expression of Eq. (21a) for the time evolution of the generating function, 
the descent equations of Eq. (7b) can be verified; this calculation is carried out in Appendix 
B. 

5. Jump process Schrodinger equation 

As our next density tensor application we consider the jump process (piecewise de- 
terministic process, or PDP) Schrodinger equation, given by 

dlV^) = Adt\^l^) + EkdNklip) , (25a) 

where a sum over k is understood, with A and the general (non-self-adjoint) operators, 
and with the dNk independent discrete random variables obeying 

dNjdNk = SjkdNk , dNjdt = . (256) 
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Straightforward calculation shows that this process preserves the norm of and the pure 
state condition = p = provided that A and B obey the restrictions 

(25c) 

{Bk + Bl + BlBk)=Q , 
with no summation over k on the second line, which must hold individually for each value 

of k. Corresponding to Eq. (25a), the density matrix obeys the evolution equation 

dp ^(Ap + pA^)dt + QkdNk , 

(25d) 

Qk^BkP + pBl + BkpBl , 
with a sum over k understood in the dN^ term on the first line, but no sum over k understood 
in the second line. 

Let now denote an expectation conditioned on the current value of the wave 

function being and E[...] be the expectation value over the entire history of the jump 
process (which leads to an ensemble of different current values of the wave function). We 
wish to find restrictions on A, Bk, and on 

E\^)[dNk] = Vkdt , (26a) 

such that the expectation of dp takes the Lindblad form of Eq. (16b), that is, 

E[dp\ ^dtCp 

1 ^ (266) 
£p = - ^[//TOT, p\ + Cfep4 - i^Wk, P} ■ 

Making the Ansatz 

= , (27a) 

with Kk constants (this Ansatz includes both the standard quantum jump equation {K^ = 0), 
and the orthogonal jump equation {K^ — (ck)), as special cases; see Schack and Brun [10] 
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for a concise review), some calculation shows that the conditions of Eqs. (25c) and (26b) are 

satisfied if we choose 

Vk ^{{ck - {ck - Kk)) , 

1 (27^) 
A = - iHTOT - ^clck + -{clck) + CkKl - -{{ck)Kl + {ck)*Kk) . 

Let us now define the order n density tensor for the jump models by 

n 

p^-^^E[l[pe] , (28a) 
e=i 

where we use the condensed notation of Eq. (5e). For the differential of this, we find 

n n 

idpi+ ^ {pi---Pn)emdpidpm 

(286) 

+ X] {pi...pn)empdp£dpmdpp + ... + dpidp2dp3...dpn-ldpn\ , 

e<m<p=l 

where all powers of dp must be retained because dN^ = dN^. Using the conditional proba- 
bility formula p( I ■0) r\dNk) — p{dNk\ \'4>))p{\'4>)) , we get the conditional expectation formula, 
valid for an arbitrary function F of the state IV'), 

E[F{\'il;))dNk] = E[F{\'il;))E\^^[dNk]] = E[F{\'il;))vk] . (29a) 

Using this equation to evaluate the higher order terms in Eq. (28b), together with Eq. (26b) 
for the leading term, we get 

n 

dp^"^^ =dtE[J2{Pi-Pn)i{jOp)e 

n 

+ XI {P^-"Pn)lmVk{Qk)l{Qk)m 



e<m=i 

n 

+ X iPl---Pn)impVk{Qk)l{Qk)m{Qk)p + ■■■ + Vk{Qk)l{Qk)2{Qk)3---{Qk)n-l{Qk)n\ , 
e<m<p=l 

(296) 

with a sum over k in each term containing Vk. 
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Writing the corresponding generating function in compact notation as 

G[a\ = E[e''-P] , (30a) 

the evolution equation for G is given , with the k sum now indicated exphcitly, by 

dG[a] ^E[e''-'^Pe''-P] - Ele"-"] 

U P- (30^) 

p=2 k P' 

Prom Eq. (30b), and the identities ( which follow, after some algebra, from Eqs. (16b), (25d), 
(27a), and (27b) ) 

{p, ^p} =^p - Yl , 

k (30c) 

{p,Qk} =Qk-Ql , 

one can prove that Eq. (30b) obeys the descent equations, as shown in Appendix C. 

6. The density tensor for quantum noise and its kinematical properties 

Let us now consider a closed quantum system, consisting of a system S interacting 
with an environment £. In such a situation, one does not have a classical probability distri- 
bution Wa and fluctuations associated with this probability distribution. Instead, one deals 
with the system plus environment as the only pure state that is given, with the fluctuations 
that are averaged over in deriving the master equation coming from quantum fluctuations 
associated with the system-environment interaction. Weighted averages of the sort that we 
have used in our definition of Eq. (4a) appear only when the total state is a mixture of pure 
states, such as a thermal state, but in this case, important system quantum fluctuations still 
occur in each pure state component of this mixture. In order to describe this more general 
situation, we shall have to generalize our definition of a density tensor hierarchy. 
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To achieve this, we initially suppose the overall system plus environment to have the 
pure state density matrix p. We denote the system basis states by \i) , as well as and 

denote the environment basis states by |ea), a = 1, 2, A general density matrix element 

has the form (eii|p|e2j), and the standard reduced density matrix, with the environment 
traced out, is defined by 

p« = (IV,p),,. = J](e^|p|ei) . (31) 

e 

In order to recapture fluctuations that are averaged over in the trace in Eq. (31), we define 
the density tensor p^") by 

(32a) 

'^^£PiljlPi2j2---Pinjn 

Here we have defined pi^j^ as the matrix, labeled by the system state labels i^, ji, acting on 
the environment Hilbert space He according to 

{Pide)eie2 ^ {ei\Pide\e2) = {eiie\p\e2je) ■ (326) 

The density tensor p*^") is again an operator on a tensor product of system Hilbert spaces 
YYe=i^S;i- Thus, in a condensed notation analogous to that of Eq. (5e), we can also write 
Eq. (32a) as 

= IV£pip2....Pn , (32c) 

where p^ is an operator acting on (8> 'Hs-g. 

We have avoided using a product notation YYe=i Eq. (32c) because the factors pi^j^ 
in Eq. (32a) and p£ in Eq. (32c) are different operators on the environment for each i and 
thus do not commute. Hence the density tensor is not symmetric under permutation of its 
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pair indices ieji, but it is symmetric under cyclical permutation of the indices, as a result 
of the cyclic symmetry of the trace. For n = 2, cyclic symmetry is equivalent to symmetry 
under pair index interchange, and for n = 3, using the identity 



cyclic symmetry is equivalent to the statement that the density tensor p(^) can be written 
as the sum of two tensors p^^^ = p^^^^ + p^^'^\ with p^^"^) completely symmetric, and p*^^^) 
completely antisymmetric, under pair index interchange. Also because the density tensor 
is not totally symmetric in its pair indices, we cannot introduce a generating function by 
imitating Eq. (5d) 

Similarly, because of factor non-commutativity, the density tensor satisfies only a 
subset of the descent equations of Eqs. (6a), (6b), and (7b). Contraction with (5^^^^ does not 
lead to a descent condition, since Si^j^pi^j^ is not unity, but rather Trsp, the reduced density 
matrix that acts on the environment when the system is traced out. Contraction of a general 
jk with a general ii for k ^ i gives nothing useful, since in general non-commuting factors 
stand between p^ and p^. However, when a column index is contracted with the adjacent 
row index ik+i, the two density matrices to which they are attached are linked to form the 
product p^ = p, and so we get the descent relation of Eq. (6b), and others related to it by 
cyclic permutation symmetry. 



As noted before, even when p^ ^ p, the descent relation corresponding to Eq. (34) is still 
useful for checking calculations. Since we cannot define a generating function as in Eq. (5d), 
in the quantum noise case we do not have analogs of the descent equations in the form of 
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TrABC = Tr^([^, B]C + {A, B}C) 



(33) 



^jli2Phji,i2j2,-,injn 




'n 



(34) 



Eq. (7b); when verifying the descent equations in the various cases considered below, we will 
work directly from Eq. (34). 

We will also consider a more general definition of the density tensor, corresponding 
to the case in which the system plus environment is in a mixed state composed of pure states 
Pa with weights Wa- Typically, a refers to an eigenvalue of a conserved quantum number of 
the total system, such as the energy; when the environment is considered in the independent 
particle approximation, with the system back reaction on the environment neglected, a then 
can refer to the energies and momenta of each environmental particle. In this case we define 
the density tensor by 

Piiji,i2j2,-,injn ~ ^^'^°'Pa;hji,i2j2,--,injn ' (SSo) 
a 

with 

Pa-}iji,i2j2,-,injn ~ Pa-Mji Poi;i2j2 ■ ■ ■ Pa;injn ■ (356) 

This definition gives information about both the quantum noise or fluctuations contained 
within each p^, and the classical noise or fluctuations associated with the probability distri- 
bution Wa- Note that in the mixed state case one could also define a density tensor that is 
a direct analog of the classical noise definition of Sec. 2, by 



njn = n ^SPc'-Mh (36) 



(n);CL 

Hjl,i2j2,---,i _ 

=1 



which would give information only about the classical noise fluctuations associated with the 
probability distribution Wa- In the examples computed in the following sections, where a 
weak coupling approximation is made, the definition of Eq. (36) typically contains no more 
information than could be gotten from a product of n reduced density matrix factors, each 
of the form Y^a'^o^T^sPc'-Mje- 
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As already noted, the density tensor p^'^^ is not measurable by any operation on 
the system Hilbert space. Its construction requires knowledge of the full system plus envi- 
ronment density matrix, which is not experimentally accessible for complex environments. 
Nonetheless is computable in any theory of the system-environment interaction, and 
we believe it to be of conceptual and theoretical interest, even if not of direct empirical 
relevance. 

We close out this section by noting that in the quantum noise case, there is no analog 
of Eq.(8a), which relates the positive semidefinite variations Vari_2 to the density tensor p*^^) 
in the classical noise case. The closest analog we find to the fluctuation formulas of Eq. (8a) 
involves the n — ?> density tensor. The reason for this is that whereas £'[1] = 1, the trace 
over the environment of unity is the dimension of the environmental Hilbert space; to get a 
unit trace over the environment we must include a factor of ps = Tr^p, the reduced density 
matrix for the environment. This pushes up the order of the density tensor involved from 2 
to 3. Specifically, let As be an operator acting on Tis, but which acts as the unit operator on 
Hs- In place of the expectations used in the classical noise discussion of Eqs. (la) through 
(2d), in the quantum noise case of system plus environment we consider the expression 
= TrspAs, which is an operator on the environmental Hilbert space. The trace of this 
operator over the environment is TrgAg = TrgTvspAs = Trs{TiC£p)As = Tr^p^^M^, giving 
the expectation of the operator Ag when the environment is not observed. On the other 
hand, the expectation of this operator formed from the environmental reduced density matrix 
is TrgpgAg. The mean squared fluctuation of this operator over the environment is positive 
semidefinite, and is given by 

TrepsiAe - TrepeAef = TrepsAl - {TrepsAef , (37a) 
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where we have used the fact that Trspg — Trp = 1. Reexpressing Eq. (37a) entirely in terms 

of the pure state density matrix p, we have 

TrsTrsp{TrspAsf - (TrsTrspTrspAs)^ 

(376) 

—^jih ^Sj2i2 ^Sjais Phji ,1232,1333 \'^hn ^83212 Piiji ,1232 ) ' 
where we have used the fact that the right hand side of Eq. (37b) involves only the symmetric 

part of the order 3 density tensor. Thus, as noted above, where a n = 2 density tensor appears 

in Eq. (8a), a n = 3 density tensor appears in Eq. (37b), and where a n — 1 density tensor 

appears in Eq. (8a), a n — 2 density tensor appears in Eq. (37b). 

7. Collisional Brownian Motion 

As our first application of Eqs. (32a-c) and Eqs. (35a,b), we consider the collisional 
Brownian motion of a massive Brownian particle immersed in a bath of scattering particles. 
We work in the approximation of neglecting recoil of the Brownian particle, and of treating 
the bath as a collection of free particles of mass m. We consider the pure state density matrix 
corresponding to definite momenta {ki} of the bath particles, calculate the corresponding 
order n density tensor defined by Eqs. (32a-c), and then average over the thermal distribution 
of the bath particles as in Eqs. (35a,b). Thus the initial density matrix for the total system, 
corresponding to the factor p^ in Eq. (33d), is 

pr^ = PiP£ , (38a) 

with the initial density matrix of the Brownian particle, characterized by its coordinate 
matrix elements {R£\p£\R^), and with ps the product density matrix for the bath particles, 

P£-Yl\h){h\ . (386) 

i 
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Since the bath particle scatterings are all independent, we focus on the effect of 

— * 

the scattering of a single bath particle, of initial momentum k, on the Brownian particle, 
which we take to be in a superposition of position eigenstates. Thus the initial state of the 
Brownian particle and the bath particle that we are considering is 



|/) = ^c^|il)|^) , (39a) 



R 

corresponding to an initial state density matrix 

Pi-\I){I\ 

= EE^«4'I^)I^)(^K^'I ■ ^^^^^ 

R R' 

The corresponding Brownian particle matrix element of pr, which is still an operator on the 
bath particle state, takes the form 

= , (39c) 

with 

p{R,R')^c^c% . (390?) 

Asymptotically, the effect of the scattering is to replace the initial state |/) by |F) = 
S\I), with S the scattering matrix. Substituting Eq. (39a), and using translation invariance 
to relate the scattering matrix S with the Brownian particle at a general coordinate, to the 
scattering matrix 5*0 with the Brownian particle at the origin, we get [11] 

1^) =S\I) = Y.c^S\R)\k) 

^ . . . . (40a) 
= J2cii\R)e~'''°'"^Soe'''°^-^\k) , 



R 
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with ^op the momentum operator for the bath particle. The corresponding final density 

matrix is then 

PF=\F){F\ 

= ^ ^ c^4, 1 R ) e-^fcoP-^^^gifeoP-^ I ^ ^ ^-ikoP-R' sle'^oP-R' ^ ^^^^^ 

R R' 

and the Brownian particle matrix element oi pp, which is again an operator acting on the 
bath particle, is 

{R\Pf\R') = p{R,Ry-'''°^-^Soe'''°^-^\k){^^^^ . (40c) 

Substituting this expression into Eq. (33d), we get 

in) 
n 

= TT p{Ri, ^^)(A?|e~'^°p-^^^(|e*^op-^«e-*^op-^^+i5'oe*^op-^^+i |A!) , 

(40d) 

n 

= JJp(^^, ^;)(^|S'(le^^°P-(^^^^+i)5o|^)e'^'('^^+i-'^^^ , 

1=1 

with i?n+i = Ri. The matrix element appearing in the final line of Eq. (40d) is one that 
is familiar from the standard calculation of the reduced density matrix (that is, P^^^,) for 
coUisional decoherence [12]. Writing 

(^|5o^e^^op-(^^^^+i)5'o|A?)e^^-(^^+^-^^) = 1 + /(^^+i - Rl^) , (41a) 

with / proportional to the square of the scattering amplitude, the product of matrix elements 
in Eq. (40d) can be written, to second order accuracy in the scattering amplitude, as 

n n 

n[i+/(^m-^;)]^i+E/(^^+i-^^) ■ (416) 

t=i e=i 
We also note that Eq. (39c), when substituted into Eq. (32c), implies that the value of p^^^ 

before the scattering is 

e=i 
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Thus when the approximation of Eq. (41b) is substituted into Eq. (40d), we get 



An) ^ _ _ 

€=1 

At this point our work is essentially finished, since the remaining steps are identical 



to the standard calculation [11,12,13] proceeding from the n = 1 case of Eq. (41d), and the 
structure of Eq. (4 Id) makes it clear how to generahze the standard result for p^^^ to the 
case of general p^"'\ In brief, the standard procedure is to multiply the right hand side of 
Eq. (41d) by the number of scattering particles, which combines with a normalizing factor of 
the inverse volume to give an overall factor of N, the scattering particle density. The effect 
of the thermal distribution ij,{k) of momenta k is taken into account by including an integral 
J dkii{k), in accordance with the mixed state procedure of Eq. (35a). Finally, expressing 
the S matrix in terms of the scattering amplitude f{k',k), and noting that the squared 
delta function for energy conservation gives an overall factor of the elapsed time, Eq. (41d) 
becomes, in the limit of small elapsed time, a formula for the time derivative of p^^\ For the 
n — 1 case, the standard answer obtained this way is 



^^-FiR-R')p('\t)^^, , (42a) 



with 



F{R) J dkp{k)^^ J dn{l - e'(^-"l^l>-^))|/(n|^|,^)|' , (426) 

where n is a unit vector which gives the direction of the scattered particle momentum 
k' — h\k\. To compare Eq. (42b) with the n — 1 case of Eq. (41d), we replace R by 
R2 = Ri and R' by R[ . Then we see that the generahzation to n > 1 is given by 

X" " " ^ "E ^(^^+1 - ^^)]P^"Ht)«,^;,...,^„^,. (42c) 

£=1 
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This is our final result for coUisional Brownian motion, giving the evolution equation 
obeyed by the order n density tensor . We see that it has the generic symmetries expected in 
the quantum noise case: although not totally symmetric in its pair indices, p^"^ is symmetric 
under cyclic permutation of these indices. As additional checks, we see that for n = 2 the 
factor involving F is 

F{R2 -R[) + F{Ri - 4) , (43a) 

which is symmetric under the interchange 1 <-> 2, while for n = 3 we have 
F{R2 - R[) + F{Rs - 4) + F{Ri - 4) ^ F^ + F^ , 

F^ =^[F(i?2 - R'l) + F{R3 - 4) + F{R, - 4) 

+F{R,-R:,) + F{Rs-R[) + F{R2-R'^)] , (436) 

F^ =^[F(4 - R[) + FiRs - 4) + ^(^1 - ^3) 

-F{Ri - 4) - F{Rs - R[) - F{R2 - 4)] , 
with F^ symmetric, and F"^ antisymmetric, under any of the pair interchanges 1 <-> 2, or 

— * — * 

1 -^-^ 3, or 2 <-> 3. Checking the descent equations is easy. Setting R'^ — R2, the term 
F(i?2 — R'l) in Eq- (42c) vanishes, so that on integrating over R[ one is left on the right hand 
side with a sum F{Ri — R'^) + F(i?3 — i?2) + ... that does not involve times 

/ , (43c) 

and so the descent equation for p*^"' (t) then implies the descent equation for its time deriva- 
tive. 
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8. The weeik coupling Born-Mcirkov approximation and the quantum optical 

master equation for the density tensor 

We turn next to the density tensor extension of the standard weak couphng Born- 
Markov approximation, that is used to give a master equation for the reduced density matrix 
p(^) for a system S interacting with an environment £. We assume a total system plus envi- 
ronment Hamiltonian i^TOT = Hs + Hs + H, with Hs and Hs respectively the environment 
and system Hamiltonians, and with H the system-environment interaction Hamiltonian. 
(We omit the customary subscript / on the interaction Hamiltonian to avoid a proliferation 
of subscripts.) We shall work in this section in interaction picture, in which the operators 
carry the time dependence associated with Hs and Hs- Thus the interaction Hamiltonian 
carries a time dependence H{t), and the density matrix obeys the equation of motion 

dpjt) 
dt 

which can be integrated to give 



^-t[H{t),p{t)] (44a) 



p{t) = p{0)-t f ds[H{s),p{s)] . (446) 

Substituting Eq. (44b) back into Eq. (44a) gives the additional evolution equation 

^^-i[H{t),pm- iys[H{t),[H{s),p{s)]] . (44c) 

One then notes that up to an error of order H^, the time argument of the factor p{s) in the 
double commutator term is irrelevant, so this factor can be approximated as p{t), giving 



dp{t) 



= -i[H{t),p{0)] - f ds[H{t), [H{s),p{t)]] , (44d) 

^0 



dt 

which is used as the starting point for the standard master equation derivation. 
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Our first step is to derive a suitable extension of Eq. (44d) for the product pip2---A 
that appears in Eq. (32c). By the chain rule, we have 



d{plp2--Pn) dpi dpi dpn s 



For each undifferentiated factor on the right of Eq. (45a) we substitute Eq. (44b), and for 
each time derivative factor we substitute Eq. (44c), with appropriate subscripts added. Let 
us now organize the terms obtained this way according to the number of factors of H that 
appear. Since Eq. (44c) contains at least one factor of H, there are no terms in Eq. (45a) 
with no factors of H. The general term in Eq. (45a) with one factor of H comes from the 
term in Eq. (44c) with one factor of H, multiphed by the product of the terms from Eq. (44b) 
with no factors of H, giving 



-i([//l(t),Pl(0)]p2(0)...p„(0)+pi(0)[i/2(i),P2(0)]...Pn(0) + ... + Pl(0)p2(0)...[//n(i),p„(0)]) . 

(456) 

The terms in Eq. (45a) with two factors of H are of two types: (1) the quadratic term in 
H on the right of Eq. (44c) times factors of p(0), and (2) the linear term in H on the right 
of Eq. (44c), multiplied by one factor of the linear term on the right of Eq. (44b), times 
factors of p(0). We now note that up to an error of order H^, in terms that already contain 
two factors of H we can replace all factors p(0) or p{s) by the corresponding p{t), since the 
differences p{t) — p{s) and p{t) — p(0) are all of order H. Collecting everything, we get the 
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following formula, which gives the needed extension of Eq. (44d) , 



dt 



-Y,Pi{t)-Pe-i{t) / ds[He{t),[He{s),pe{t)]]pe+i{t)...pn{t) 

~ y~l'f'^l(^)---'^^-l(^)[-^^(^)''^^(^)]'^^+l(^)--'^"»-l(^) / (^4Hm{s), Pm{t)]Pm+l{t)---Pn{t) 

Jo 



+ 



Pl{t)...pe-i{t) f ds[He{s),pe{t)]pe+i{t)...prn-l{t)[H,,{t),Prn{t)]Prn+^^^^ 

Jo 

(45c) 

Taking the overall Trg of this expression then gives a formula for the time evolution of p*^"^ (t) 

as defined by Eq.(32c). 

We now make two standard assumptions. First of all, we assume at that at the 

initial time t = 0, the density matrix factorizes so that p(0) = peps, with ps and ps 

respectively density matrices for the environment and the system which commute with one 

another, and with ps a pure state density matrix obeying p| = p£. Secondly, we assume 

that {H)£ = TtsPsH = 0, that is, we take the interaction Hamiltonian to have a vanishing 

expectation in the initial environmental state. As a result of these two assumptions, the 

environmental trace of the first term on the right hand side of Eq. (45c) vanishes, since 

Tr£p,{0)...pt.i{0)[He{t),pe{0)]pe+i{0)...pn{0) 

(46a) 

^PSl---PSi-l[i^r£PsHi{t)),psi\psi+l---PSn = . 

The remaining terms in Eq. (45c) all have two factors of H. Since p{t) and p(0) differ by 
one power of H, in these terms, up to an error of order H^, we can replace all factors p{t) 
by the factorized approximation 

p{t) ~ p(0) = Peps = PsTrepiO) ~ peTrepH) = pep^'^ (t) . (466) 

With these simplifications, and remembering that system operator factors p^^"* with different 
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index values £ act on different Hilbert spaces TCs-/ and so commute, Eq. (45c) becomes an 
extended version of the Redfield equation, 



X {[Heit),pf\t)pe]pr'-'[Hmis),p';!,\t)ps] + ms), pPit)pe]pr'-'[Hmit), p';l,\t)pe]} 



This is converted to the Born-Markov equation by setting s ^ t — s, and then extending the 
upper hmit of the s integration from i to oo, giving 



X {[H,{t),pf\t)pe]pr'-'[Hm{t - s),p^^\t)ps] + mt - slpf\t)ps]pr'-'[H^it),p^(t)pe]} 

(46rf) 



We now note that Eq. (46d) can be further simphfied, by taking account of the fact 
that whenever an H factor is sandwiched between factors of ps it vanishes, since psHpg — 
Ps{H)e = 0. This ehminates all terms in the sum over £, m that are not adjacent in a cyclic 
sense, i.e., that do not either have m — I -\- 1, I — 1, n — 1, or £ — l,m — n. The latter, 
by use of the cyclic properties of the trace, can be rearranged to give the £ = n term of the 
former set. We thus get a simplified set of Born-Markov equations. For n = 1, we get the 
usual starting point for the Born-Markov master equation derivation. 




dp^-\t)/dt ^-J2{p?{t)...p^^\t))eTTePr' / ds[H,{t), ms),pf\t)ps]] 




(46c) 




n 



dp^^^ {t)/dt = -Trg / ds[H(t)H(t - s)p^^^ {t)p£ + p(i) {t)p£H{t - s)H{t) 




(47a) 



-H(t)p^'\t)p£H(t -s)- H(t - s)p^'\t)psH(t)] 
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and for n > 2, with the subscript n + 1 identified with 1, 

POO 

{t) jdt = -Tl£P£ / dsJ2 

Jo 

X {{p^^\t)...p^:\tMH,{t)H,{t - s)pf\t) + pf\t)H,it - smt)] 

- {p?it)...pl^Ht))u+i[p?m,it)He+^{t - s)pf^,{t) + pf\t)H,{t - s)He+^it)pflm . 

(476) 

At this point it is useful to check (and we have done so) that the descent equations are 
satisfied by Eqs. (47a) and (47b). 

The remainder of the derivation foUows closely the standard master equation deriva- 
tion, in the rotating wave approximation, that proceeds from Eq. (47a), so we will only give 
a sketch. For further details, and in particular a discussion of the physical justification for 
the approximations involved, see Sec. 3.3 of ref [1] and also ref [13]. One assumes that H£(t) 
has the form 

a CO 

with acting only in the system Hilbert space Hs;e and with acting only in the 
environment Hilbert space He, and with the Hermiticity properties AI^{uj) = A^q,(— cu) and 
Bl^{t) — Ba{t). Since Eqs. (47a,b) are quadratic in H, one uses Eq. (48a) twice; for each 
Hk{t — s) (regardless of the value of the index k) one writes 

H,{t-s) = J2 e-'^^'-'^A,fs{uj)Bf,{t - s) , (486) 

l3u) 

and for each Hk{t) (again regardless of the value of k) one writes 

au>' 

The rotating wave approximation then consists of neglecting terms in the double sum with 
a;' 7^ a;, so that only the diagonal terms cu' — cu are left. Prom the trace over the environment, 
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and the integral over s, one gets correlators of the form 

poo 

/ dse'^'{Bl{t)Bf,{t-s))£ = r^p{u;) , 

-^'oo (49a) 
/ dse'^^{Bp{t-s)Bi{t))£ = T^p{-ujr , 

where in the second line we have used the definition of the first line and the adjointness 
properties of the integrand. It is also customary to decompose the reservoir correlation 
function Tap into self-adjoint and anti-self-adjoint parts, according to 

ra/3(c^) = ^7a/3M+«'5a;3M ■ (496) 

Proceeding in this fashion, after some algebra one gets the final result, which can be 
written as an equation for all n > 1 by including a 5ni to take account of the special nature 
of the n = 1 equation, 

n 



n 



+ EE^"/^H \ iP?it)...p(:\t))e S^,Ae,{u;)p^l\t)Al{u;) - i{AL(a;)A,,(a;), p«(t)} 
+ (p«(t)...p«(t))«-,ipW(t)4,(a;)A,+i,(a;)pa(t)} . 

(50a) 

Despite the fact the the n = 1 and n > 2 density tensors have a different structure, the 
descent equations are satisfied by Eq. (50a), as verified in Appendix D. 

Finally, we note that Eq. (50a) is readily converted to the quantum optical master 
equation and its density tensor generalizations, by taking a to be a three-vector index, 
so that Aa becomes A, which is related to the dipole operator by Eq. (3.182) of ref [1]. 
Also, one takes Sap{(^) = dapSiuj), with S{uj) given by Eq. (3.205) of ref [1], and "^apij^) = 
(4a;^/3)[l -|- N{u;)]Sap, with N{u;) — l/{e^^ — 1) the photon number operator. One gets in 
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this way the density tensor generahzation of the quantum optical master equation, 

n 

£=1 oj 
n . 

+ EE(4^V3)[1 + A^H] {P?{t)...p^:\t)) 



(506) 

which is our final result of this section. 

9. The Caldeira-Leggett model master equation for the density tensor 

The Caldeira-Leggett model [14] describes the damping of the one-dimensional mo- 
tion of a Brownian particle of mass m, moving in a potential V{x), and interacting with 
an environment consisting of harmonic oscillators with masses rrio and frequencies cuo, and 
annihilation operator bo- The interaction Hamiltonian is assumed to be a linear coupling 
H — —xB, with 

B = y^KpXo = y^^Ko{bo + bl)/{2moU}o)^ (51a) 

o o 

a weighted sum of the harmonic oscillator coordinates. A counter-term formally of order H^, 

H^^x^Y ^^=x''C , (516) 

o 

is included in the calculation, so that the total Hamiltonian is 

Htot^Hs + Hs + H + Hc , (52a) 



with and Hs respectively the oscillator and particle Hamiltonians, 

Hs^j2^o{bibo+h , 



o 

-.2 



(526) 
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Our aim will be to get a description of the effect on the particle motion of the couplings 
to the oscillator environment, in the high temperature limit. Our derivation of the density 
tensor generalization of the high temperature master equation closely follows that of Sec. 
3.6 of ref [1], to which the reader is referred for a discussion of the physical motivation of 
the approximations involved. 

Since the environmental expectation of the interaction Hamiltonian H vanishes, we 
can proceed directly from the simphfied Born-Markov equation of Eqs. (47a) and (47b). The 
first step is to transform the density matrix p{t) back to Schrodinger picture; it is easy to 
see that the effect of this is to replace H{t) hj H = H{0), to replace H{t — s) by H{—s) 
(with H{—s) still in the interaction picture), and to change d/dt to D/dt, defined by 

n 

Dp^^\t)/dt = dp('^)(i)/dt + i^Trfpi(i)...p,_i(t)[p,V(2m) +y(x,),p,(i)]p,+i(i)...p„(i) . 

(53a) 

It is also necessary to explicitly include commutators arising from the counter term, which is 

easy since this term is treated as being already quadratic in H. For the analog of Eq. (47a) 

for the special case n = 1, we find 

Dp^^^ {t)/dt = -i[H,, (t)] - Tte / ds[HH{-s)p^'^ {t)p£ + p^'^ {t)psH{-s)H 

•^0 (536) 

- {t)peH{-s) - H{-s)p^'^ {t)peH] , 
and for the analog of Eq. (47b) for n > 2, we have 

n 

dp^-\t)/dt = -lY,{P^\t)...p^nKmHci,pf\t)] 

/•oo ^ 

-T^sPs / dsY, 
Jo 

-(p«(t)...p«(0)«+i[p?^(t)//.//m(-«)p2i(^) + pf\t)He{-s)He^JllM} ■ 

(53c) 
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We next note that 

Hi = -xe{0)B{0) , Hi{-s) = -Xi{-s)B{-s) , (54a) 

where, using the assumption that the system evolution is slow compared to the oscillator 
time scale, we approximate Xi{—s) by its free particle dynamics, 

Xi(-s) c:txi-—s . (546) 
m 

Since the right hand sides of Eqs. (53b, c) are quadratic in H, the operator B giving the 
coupling to the oscillators appears, after the environmental trace is taken, only through the 
correlators 

D{s)^i{[B{0\B{-s)])e , 

(55a) 

D^{s)={{B{0),B{-s)})£ , 

so that we have 

{B{0)B{-s))e=]:[D,{s)-zD{s)] , 

' (556) 
{B{^s)B{0))e=-[D,{s)+tD{s)] . 

These correlators appear in the following integrals, which are evaluated or approximated in 

Sec. 3.6.2 of ref [1], 

POO 

/ dsD{s) =2C , 
Jo 

POO 

/ dsDi{s) =Am'ykBT , 

(55c) 

dssD{s) —2mj , 



Jo 



'0 
poo 

/ dssDi{s) ^Am'jkBT/Q ~ , 
Jo 

with C the constant defined by the counter term of Eq. (51b), with 7 a constant determined 
by the harmonic oscillator spectral density, with ks and T respectively the Boltzmann con- 
stant and environment temperature, and with Q a frequency cutoff. For a spectral density 
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J(a;) with a Lorentz-Drude cutoff function, one has 

o 

This completes the specification of the calculation; the rest is just the algebra of 
assembling all the pieces, and so we pass directly to the result. For n = 1, we get the 
Caldeira-Leggett master equation, 

Dp^^\t)/dt = -i-f[x, {p,p^^\t)}] - 2m-fkBT[x, [x,p^^\t)]\ . (56a) 

For the density tensors with n > 2, we correspondingly get 

Dp^-^/dt=j2iP^^\t)...p^:\t))e 
e=i 

n 

+ ^(P? (^)--Pn ^ {'t))ie+i [4:m'ykBTpf^ {t)xiXe+ipfl^ (t) + i'ypP (t) {xipi+i - PiXt+i)p\ 
e=i 

(566) 

We also note that the term proportional to 7 on the first line of Eq. (56b) can be written in 
the alternative form, 

n{pf\t)pexe - xepepf\t)) = -fpf\t) + ^7[p?V)> {xe,Pe}] ■ (56c) 

Equations (56a) and (56b) are our final results for the Caldeira-Leggett model. As was 
the case for the master equations derived in the preceding section, despite the differences 
between the structure of the n — 1 and the n > 2 equations, the descent equations are 
satisfied, as verified in Appendix E. 

10. An application to state vector reduction 

We turn now to considerations that bridge the discussions given above in the classical 
and quantum noise cases. We begin with an analysis of two Ito stochastic Schrodinger 
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-2m-fkBT{xl, p[^\t)] + i-i[p^l\t)pixi - xepepf\t)) 



equations, 

1 



d\^l;)^--{A-{A)Ydt\^l;) + {A-{A))dWt\^l;) , (57a) 
and 

= -\A^dt\'il)) + iAdWt\il}) , (576) 

with dWt a real Brownian noise obeying dW^ — dt, and where we have dropped the Hamil- 
tonian term. These lead to the respective density matrix evolution equations 

dp^ ~[A,[A,p]]dt+[p,[p,A]]dWt , (58a) 

and 

dp^~[A,[A,p]]dt + i[A,p]dWt , (586) 

which correspond to the same Lindblad type evolution equation for the expectation E[p\, 

dE[p] ^ jCE[p]dt , /:p^~^[A,[A,p]] . (58c) 

Let us now consider the effect of the stochastic evolutions of Eqs. (58a,b,c) on the 
expectation of the variance V — Var(A) of the operator A, 

V =TrpA^ - {Tr pAf , 

E[V] ^TrE[p]A' - E[p,^j^p,^j,]Aj^,^Aj,,^ (59a) 

where in the final hne we have used the density tensor definition of Eq. (15). For the time 
evolution of E[V] we have 

dE[V]/dt =Tt{CE[p])A' - dp^l^,^^^A,,,,A,,,, 

(596) 

—Tr{CE[p\)A — 2E[pi^j^{Cp)i2j2]Aj^i^Aj2i2 — ElCi^j^^i^jr^Aj^i^Aj^i^ , 
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where we have used Eq. (58c) in the first hne and Eq. (23a) in the second hne. Since the 
cychc property of the trace implies that Tt[A, [A, p]\A = , the terms in Eq. (59b) involving 
the Lindblad C all vanish, and so the time derivative of E[V] comes entirely form the final 
term, 

dE[V]/dt = -E[Ci,j,,i,j,]Aj,i,Aj,i, , (59c) 

and thus is determined by the evolution equation for the second order density tensor. This 
is why the state vector evolutions of Eqs. (57a) and (57b), or equivalently the density matrix 
evolutions of Eqs. (58a) and (58b), lead to very different results for the evolution of the 
variance of the operator A. The tensor Ci^j^^i^j^ corresponding to Eqs. (57b) and (58b) is 
given in Eq. (24b), and since the cyclic property of the trace implies that Tr[A, p\A — 0, one 
has dE\y]/dt = for this evolution. On the other hand, the tensor Ci^j^^i^j^ corresponding 
to Eqs. (57a) and (58a) is given in Eq. (24a), and through Eq. (59c) implies that 

dE[V]/dt = -E[{Tt[p, [p, AUY] = -E[{Tt{[p, A]m , (59d) 

which is negative definite. Starting from Eq. (59d), some simple inequalities imply that the 
stochastic evolution of Eqs. (57a) and (58a) drives the variance of A to zero as t — > oo, and 
hence reduces the state vector to an eigenstate of ^4, as discussed in detail in refs [15]. 

Let us now consider a quantum system S, consisting of a microscopic system coupled 
to a macroscopic measuring apparatus, interacting with a quantum environment S, with the 
totality forming a closed system. A general result [16], using just the linearity of quantum 
mechanics, shows that state vector reduction cannot occur in this case. To understand this 
result through an analysis similar to that just given for Eqs. (57a,b), let us consider the 
behavior of the variance of a system operator A which is a good "pointer observable" . By 
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definition, a system operator commutes with the environment Hamiltonian Hg, and since 
the system in this case includes the apparatus and so is macroscopic, the pointer observable 
also obeys [17] [A, H] — 0, with H the system-environment interaction Hamiltonian. Let us 
now write the density matrix evolution in Schrodinger picture, 

dp/dt - -i[HTOT: p] = -i[Hs + H£ + H,p] , (60a) 

with Hs the system Hamiltonian. We consider the system evolution after a brief interaction 
has entangled the apparatus states with the microscopic subsystem quantum states that are 
to be distinguished by the pointer reading. For the time evolution of the variance of the 
pointer observable A, we have 

dV/dt = Tr{dp/dt)A^ - 2{TrpA){Tr{dp/dt)A) , (606) 

which substituting Eq. (60a), and using the cyclic property of the trace and the fact that A 
commutes with both and H, simplifies to 

dV/dt = iTrp[Hs, A^] - 2i{TTpA){TTp[Hs, A]) . (60c) 

This can be further simplified by using the definition of the reduced density matrix p^^^ — 
Trgp, together with the fact that the commutators in Eq. (60c) involve only system operators, 
giving 

dV/dt = iTrsp^'\Hs,A''] - 2i{Trsp^'^ A){Trsp^'^[Hs, A]) . {6Qd) 

We see that, unlike the Ito equation case discussed above, the time derivative of V here is 
determined by p^^^ , rather than by p^^^ . 

Let us now take the pointer observable to be a pointer center of mass coordinate 
A — X, in which case, once the entanglement of the pointer with the microsystem being 
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measured has been established, the relevant part of the system Hamiltonian Hs is P^/ (2M), 
with P the total momentum operator for the pointer of macroscopic mass M. Evaluating 

the commutators, and writing Tr^p^^^O = (C), we see that 

dV/dt ={1/M){{P,X}) - (2/M)(X)(P) 

(61a) 

^{1/M){{X-{X),P-{P)}) . 
By the Schwartz inequality, the right hand side of Eq. (61a) is bounded by 

{2/M){{X - {X)f)-2{{P - (P))')^ = {2/M)AXAP . (616) 

Let us now determine the minimum value of the bound of Eq. (61b) that is compatible 
with the parameters of a feasible measurement. Since the uncertainty principle implies that 
AXAP > 1/2, we get a least upper bound on Eq. (61b) by substituting AXAP ~ 1/2. 
This shows that \dV/dt\ can be made as small as ~ 1/M, which since M is macroscopic, can 
be made essentially arbitrarily small.^ Hence the variance of the pointer variable A stays 
essentially constant, and is not forced to reduce to zero in the course of the measurement. 

We conclude, in agreement with the arguments of [16], that a quantum apparatus 

interacting with a quantum environment does not act like the stochastic equation of Eq. (57a) 

in terms of reducing the state vector. Although a quantum environment acts on a quantum 

system with a form of "noise" , our analysis of the density tensor hierarchy in the classical and 
^ Restoring factors of Planck's constant, \dV/dt\ can be as small as k/M, for which 

the reduction time dt is at least of order MdV/h. For M ~ lO^^mproton and dV ~ (Icm)^, 

this gives dt ~ M(lcm)^/fi ~ lO^'^'s ~ 10^° times the age of the universe. Note that our 

argument places no restriction on the mean pointer momentum (P) that establishes the time 

needed to attain one or the other of the measurement outcomes X starting from the initial 

pointer position. 
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quantum noise cases shows that structures with different kinematical symmetries,^ different 
dynamical evolutions, and different implications for the measurement process are involved. 
As a result, the quantum noise in a closed quantum system does not mimic the action of the 
classical noise in objective reduction models, and cannot be invoked to give a resolution of the 
quantum measurement problem within the framework of unmodified quantum mechanics. 
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^ The dissimilarities between the symmetries of the classical noise and quantum noise 
hierarchies are least for the order two density tensor. In the order two case, cyclic symmetry 
is equivalent to full permutation symmetry, and so the index symmetry properties are the 
same in the classical and quantum noise cases, and as a consequence the descent equations 
in the quantum noise case correspond to the idempotence descent equations in the classical 
noise case. Only the classical noise descent equation implied by the unit trace condition has 
no precise quantum noise counterpart: in the classical case, one has 

whereas in the quantum noise case one instead has 

^nhpfih,i2j2 ^ ^iljl'^ePhjlPi2j2 = T^eP£Pi2j2 T^£Pi2j2 = Pil]^ ' 

with p£ — Tr sp the reduced density matrix of the environment with the system traced out. 
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Appendix A: Descent equations for the 
isotropic spin- 1/2 ensemble 

Let us write the generating function of Eq. (14b) as 



/(x) = sinhx2 /x2 , x = A'^ , (^41) 



Then, we find 



ddrn.r 2 



"mr 

32 



d G 1 dG 1 J. „/ 7* / Ar,\ 

Omr^ V 7;0pqgj A ■ Gmr K^^) 



Here the primes denote derivatives of / with respect to and in this notation / obeys the 
second order differential equation 

Contracting the first expression in Eq. (A2) with 8mr-i and using the tracelessness of the Pauli 
matrices, gives the first equation in Eq. (7b). Contracting the second expression in Eq. (A2) 
with drp-i and using the differential equation of Eq. (A3) together with the Pauli matrix 
identities {a'^)mq = •iSmq and crV-' = 5^^ + ie^^^a^, which implies A - dmpA - a^q = A'^6mq, gives 
the second equation in Eq. (7b). 

Appendix B: Descent equations for the 
Ito stochastic Schrodinger equation 

We wish here to verify that 

dG[a] = dtE[(amr{jO-p)mr + ^amrapqCmr,pq)e'''"'] {Bl) 
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obeys the descent equations of Eq. (7b). Since 

^mr^^P)mr ^mrC'mr,pq ^pqC'mr,pq ; (-B2o) 

we have 

^uv^ = dtE[{a„,rijCp)mr + 7;(^mrapqCmr,pq) i^P)^^ ''] = dG[a\ , (526) 

giving the first identity in Eq. (7b). Next we calculate 



{B3a) 



2 

while for the contraction of the second variation we have (with indices m, q implicit on the 



right hand side) 

d^dG[a] 



dtE [ (5i + + a^v {Tiuv + T^uv) ) e''-"] , (S36) 



with 



'S'l — 2 i.^mr,rq ~l~ Crq^mr) i 



mq 1 



{B3c) 



Tluv — 2 [i^'mr,uv ~l~ Cuy^^rnr^Prq ~l~ Pmr{Cuv,rq ~l~ C'j-g^u^,)] , 
T2UV [(-^P)?*!) ~l~ '^^rsCuv,rs\Pm,q ) 

We see immediately that auvT2uv gives all of the second line of Eq. (B3a). Prom Eqs. (16b) 
and (18) we find 

{Cp, p}mq = {J^p)mq " [(Cfe - {Ck))p{Ck - (Cfe))'^]mg " Pmq{{Ck " {Ck))\ck - {Ck))) , (54a) 

while from Eq. (21b) we have 

^{Cmr,rq + Crq,mr) = [(Cfe - {Ck))p{Ck - (Cfe))'^]mg + Pmq{{Ck - {Ck))\ck - {Ck))) ■ {B4b) 
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Hence Si + S2 — {Cp)mq, giving the Cp part of the first fine of Eq. (B3a). Finally, again 
using Eq. (21b) we find that 

and so auvTiuv gives the remainder of the first line of Eq. (B3a), completing the check of the 
descent equations. 

Appendix C: Descent equations for the 
jump Schrodinger equation 

We verify here that 



p=2 k ^ 



obeys the descent equations of Eq. (7b). Since Tr(£p) = and TrQk = (Bk+Bl+BlBk) = 0, 



we have 



S^.^^^dtE[{a.Cp + J2T.^k^-^^){Trp)e-'^] , (Clb) 



checking the first line of Eq. (7b). Next we calculate the first variation of G, 

P 

and the contracted second variation, 

d''dG[a] 



(C2a) 



pi 

p=2 k ^ 



dtE [{Si + 52 + ^3 + Si)e''-P] , {C2h) 
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with 



p=2 k 
S2 ={>Cp, p}mq , 



00 



p=2 k ^ 

\ P=2 k ^' J 

We see immediately that 5*4 gives all of the second line of Eq. (C2a). From Eq. (30c), which 
we rewrite here, 

{^P, P} ^^P - Yl ^^^fe ' 

k (C3a) 

{Qk,p}^Qk-Ql , 

we see that the Cp part of S2 and the Qk part of {Qk, p} in give the first hne of Eq. (C2a). 
To complete the verification, we must show that Si cancels against the remainder of 82 + 83, 
which is 

k p=2 k 

But separating off the p = 2 term of 5*1, and making the change of variable p — > p + 1 
in the remaining sum, we see that 81 is exactly the negative of Eq. (C3b), completing the 
argument. 

Appendix D: Descent equations for the Born-Mcirkov 

master equation 

We wish here to verify that Eq. (50a) obeys the descent equations of Eq. (34). We 
separate the verification into two parts, first checking the descent from n = 2 to n = 1, 
and then checking the descent from general n > 2 to n — 1. For the n — 2 density tensor 
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time derivative, writing out all terms in Eq. (50a) explicitly, and using the fact that since 
operators labeled with subscripts 2 and 1 act on different Hilbert spaces, the order in which 
they are written is irrelevant, we have 



(Dla) 

Contracting the column index associated with the subscript 1 with the row index 
associated with the subscript 2, and dropping the subscripts since all operators now act in 
the same Hilbert space, we get 

dp^'\t)/dt ^t[p^'\t),J2Sa^i^)AiHM^)]P^'\t) 



u>a/3 

'2' 



- E lap (^)p(^^ {t) \ {Ai (^)^M^) , p^'' m (Dib) 

+ E7a/9H[p^^H^)^iM^/?(^)p^'H^) + ^/3H(p(^H^))'^iH] 

^t[{p'^'\t)r,J2Sap{u^)Ai{u;)A,{u;)] 

which has the structure of dp^^\t)/dt and so verifies the 2 — > 1 descent. 
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To verify the n — > n — 1 descent we make some simplifications in notation. We omit 
all superscripts (1), since this leads to no ambiguities, as well as all time arguments (t) and 
all frequency arguments (cu). We also abbreviate 

{D2a) 

Our general strategy is to split the sum Yl^=i containing {p\...pn)i into 'Y^^Zl plus the I — 1 
and the £ = n terms, and to split the sum Yl]^=i containing (pi...pn)«+i into YyeZ2 P^^^ 
£ — 1, £ — n — 1, and £ — n terms. For the part of dp^'^^dt involving L^, we have 

n-l 

^ipi...Pn-l)iPni[Pl, Li] + {p2--Pn)i[Pl, Li] + {pi...pn-l)i[Pn, Ln] , {D2h) 

which on contracting the column index associated with the subscript n with the row index 
associated with the subscript 1, and relabeling all quantities that had subscript n with 
subscript 1, since they act now in the same Hilbert space, gives 

n-l 



^{p\p2...pn-l)ii[pt, Le] + P2--Pn-li{Pl[Pl: Li] + [pi, 

e=2 

n-l 

^ ^{plP2-Pn-i)ei[pi, Le] + p2-Pn-ii[pl, Li] , 



1=2 



which has the correct structure for the corresponding part of dp^'^ ^^dt, with pi replaced by 
pI. The remainder of dp'^^^^ /dt is 



n-l 



[1 ^ 11 1 

- ^{pi...pn)i-{Mi,pi] - (p2...p„)-{Mi,pi} - (pi...p„-i)-{M„,p„} 

i=2 

(n-2 
^{pi-Pn)ee+iPeAl^Ae+ippe+i+ p3...pr^piA\^A2pp2 
1=2 

+ P\--Pn-2Pn-\A^n-\(x^n^Pn + P2-- Pn— 1 Pn 
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Again, contracting the column index associated with the subscript n with the row index 
associated with the subscript 1, and relabehng all quantities that had subscript n with 
subscript 1, since they act now in the same Hilbert space, gives 

n-l ^ / I \ 

- J2ipi-Pn-l)i^{Mi,Pi} - (p2...Pn-l) ( PlMipi(*) + -{M,,pI}] 
1=2 ^ ' 

(n-2 
^{pl-Pn-i)u+ipeAl^Ae+ippi+i + {p3...pn-i)plA{^A2pp2 
1=2 

+ p2---Pn-2Pn- P2---Pn—1 

which on canceling the terms marked with (*) gives the corresponding part of dp^"'^^/dt, 
with pi replaced by pf. This completes the verification of the n — > n — 1 descent. 



Appendix E: Descent equations for the 
Caldeira-Leggett model 



We verify here that Eqs. (56a) and (56b) obey the descent equations of Eq. (34). 
As in the preceding appendix, we simplify the notation by omitting all superscripts (1) and 
all time arguments (t). We first verify the n — 2 to n — 1 descent. For the n — 2 case of 
Eq. (56b), we have 

Dp^'^^/dt =p2[-2mjkBT{xlpi + pixj) + ij{piPiXi - Xipipi)] 
+pi[-2m-fkBT{xlp2 + P2xl) + ij{p2P2X2 - X2P2P2)] 

+^m'ykBT{piXiX2P2 + P2X2X1P1) + i-f[pi{xiP2 - Pi 2^2)^2 + ^2(2^2^1 -^22^1)^1] • 

{ElA) 

Contracting the column index associated with the subscript 1 with the row index associated 
with the subscript 2, and dropping subscripts since all operators now act in the same Hilbert 
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space, we get 

Dp^'^^ jdt — > — 2m'~)kBT{x'^ + px^ p) + i'y{ppxp — xpp^) 

— 2m'ykBT{px'^ p + p^x'^) + i^{p^px — pxpp) (ElB) 

-\-A.m^kBT[px^ p + xp^x) + i^[p{xp — px)p + pp^x — xp^p] 
We see that the terms that have an operator sandwiched between two factors of p cancel, 

leaving only terms involving p^, which have the form of Eq. (56a) with p replaced by p^. 

To check the n > 2 to n — 1 descent, we split the sums that occur in the same manner 

as in Appendix D. We thus write Eq. (56b) in the form 

n-l 

Dp^'^'^/dt = ^{pi...pn)t[-'2m-ikBT{x1, pe} + i-f{pipexi - xiptp^)] 

1=1 

+P2--Pn[-2m'ykBT{xj, pi} + i'y{piPiXi - Xipipi)] 

+Pl...Pn-l[-2m'jkBT{xl, Pn} + iliPnPnXn " X^PriPn)] 
n-2 

+ ^(Pi-Pn)u+i[4^rn'ykBTpiXeXi+iPi+i + i-fPi(xePi+i - pexe+i)pi+i] ^"^^^ 
e=2 

+Ps...Pn[^m'ykBTpiXiX2P2 + hPi{xiP2 -piX2)p2] 
+Pl---Pn-2['imjkBTpn-iXn_iXnPn + ijPn-l{x„-lP„ - P„-lXn)Pn\ 
+p2---Pn-l['^m-fkBTpnXnXiPi + ijpniXnPl " PnXl)pl] • 

We now contract the column index associated with the subscript n with the row index 
associated with the subscript 1, and relabel all quantities that had subscript n with subscript 
1, since they act now in the same Hilbert space. As is readily seen by inspection of Eq. (E2), 
this gives Eq. (56b) with n replaced by n — 1 and with pi replaced by pf, together with terms 
of the wrong structure, that grouped together give (4 — 2 — 2)p2...pn-\rn'ykBT pix\pi — and 
(1 — l)p2...pn-iilPi{xipi — piXi)pi = 0, which thus vanish. This completes the verification 
of the descent equation for Eq. (56b). 
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